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Abstract 

We study the energy spectra of small three-dimensional (3D) and two- 
dimensional (2D) semiconductor quantum dots through different theoreti- 
cal approaches (single-site Hubbard and Hartree-Fock hamiltonians); in the 
smallest dots we also compare with exact results. We find that purely 2D 
models often lead to an inadequate description of the Coulomb interaction 
existing in realistic structures, as a consequence of the overestimated carrier 
localization. We show that the dimensionality of the dots has a crucial impact 
on (i) the accuracy of the predicted addition spectra; (ii) the range of validity 
of approximate theoretical schemes. When applied to realistic 3D geometries, 
the latter are found to be much more accurate than in the corresponding 2D 
cases for a large class of quantum dots; the single-site Hubbard hamiltonian 
is shown to provide a very effective and accurate scheme to describe quantum 
dot spectra, leading to good agreement with experiments. 
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I. INTRODUCTION 



Adding an electron into a semiconductor quantum dot (QD) produces a variation in the 
energy of the system that depends on single-particle quantum confinement as well as on 
the Coulomb interaction between carriers.0 Understanding such addition- energy spectrum 
is a key step towards controlling the physics of single-electron devices. At the same time, 
the addition spectra of quantum dots offer a unique probe of few-particle interactions in 
regimes that are not experimentally accessible in atomic physics. The experimental effort in 
this direction has developed very rapidly after the recent fabrication of controlled small-QD 
devices based on gated vertical heterostructuresi or self-assembled dots.! The resulting ad- 
dition spectra show a clear shell structure, corresponding to the symmetries of the confining 
potential, with a filling sequence analogous to Hund's rule in atomic physics. 

From the theoretical point of view, a general interpretation of these features was obtained 
by calculating the energy spectrum for a strictly two-dimensional (2D) quantum dot, and 
using either exact methods (for very few electrons), or approximate — usually Hartree-Fock — 
methods.^ The assumption of a purely 2D model was initially motivated by the typical disk- 
like shape of the QD potential, whose extension along z is (slightly) smaller than the lateral 
extension of the carrier ground state in the xy plane. If one adopts a separable picture for the 
QD confining potential, V = V(z) + V(x,y), the relevant (i.e. lowest) single-electron states 
can be all associated to the ground state of V(z). From the point of view of single-particle 
states the 2D assumption is therefore justified. 

In view of the three-dimensional (3D) nature of the Coulomb interaction, however, the 
2D model introduces additional approximations in the calculation of the Coulomb integrals, 
which are sensitive to the spatial extension — 2D vs. 3D — of the single-particle wave 
functions.!!! In turn, Coulomb integrals control electron-electron correlation, and influence 
the quantitative determination of addition spectra and their dependence on magnetic field. 
At the same time, the strength of Coulomb interaction is also the key parameter determining 
the accuracy and range of validity of the approximations which must be introduced for dots 
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with many electrons. 

In this paper we investigate theoretically the addition spectra of realistic QD structures, 
with special emphasis on the effects of electron-electron repulsion and their dependence on 
the geometry and dimensionality of the confining potential. In Sect. [TI], we compare different 
approximate solutions of the general Hamiltonian for N interacting electrons confined in a 
QD structure; in particular we consider the single-site Hubbard (SSH) scheme introduced 
in Ref. ||, and the standard Hartree-Fock (HF) method. 

In Sect. [Ell], we focus on the simplest CcLSG, 1.6. . cL two-electron system within a parabolic 
confining potential, and calculate the exact energy eigenvalues and pair correlation functions 
for the 2D and 3D case. As in Ref. [7], we use this prototypical system — called artificial or 
QD Helium — as a reference to evaluate the accuracy of the different approximation schemes: 
We find that both the importance of corrections beyond HF and the differences between HF 
and SSH are drastically reduced for a realistic 3D description of the dot with respect to its 
2D modelization, mainly as a consequence of the reduced Coulomb integrals. This suggests 
the reliability of a fully 3D mean field treatment of semiconductor QDs. 

Section [TV] is then devoted to the application of HF and SSH methods to 3D and 2D 
quantum dots with a larger number of electrons. We compare both methods for QD struc- 
tures of different geometries and demonstrate that SSH is an accurate and efficient scheme 
for realistic, i.e. 3D-like, dots. Finally, we discuss the implications of our results for the 
interpretation of recent experimental data vs magnetic field in QD structures and draw some 
conclusions. 

II. THEORETICAL APPROACH: EXACT FORMULATION AND 
APPROXIMATION SCHEMES 

Our aim is to describe N electrons, confined in a QD structure (with harmonic in-plane 
confining potential) and interacting via Coulomb law, possibly in the presence of an external 
magnetic field perpendicular to the plane. The general iV-particle Hamiltonian is 
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where the single particle Hamiltonian, within the effective mass approximation, is 



^o(0 = i + 7 A (^)) 2 + + y 2 ) + V( Zi ) . (2) 
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Here A is the vector potential, k and m* are the scalar dielectric costant and the effective 
electron mass in the semiconductor, ujq is the characteristic oscillator frequency of the in- 
plane confining potential, and V(z) is the confining potential along z; V(z) can be chosen 
either as a harmonic potential (V(z) = ^m*UQ Z 2 ), a square well or a zero-width infinite 
barrier to describe spherical, cylindrical or disk-shaped QD structures, respectively. Here 
Zeeman coupling between spin and magnetic field has been neglected. 

This general Hamiltonian can be written in second quantized form on the complete and 
orthonormalized basis of single particle states 

1 

Here e a are the eigenenergies of the one-particle Hamiltonian H , & aa and c aa the creation 
and destruction operators for an electron with orbital index a and spin a; V aa ^ c i.^ a i ^ a are 
the two-body matrix elements of the electron-electron interaction 

Fa*( r > s ) ^( r '> s> ) I _ n <t>iA r \ s ') Mr, s) drdr' 

SS' I I 

where (j) aa (r, s) = <f> a (r) Xo-( s ) are the single particle eigenfunctions. 

It is useful to isolate among the Coulomb matrix elements the "semi-diagonal" ones, 
namely 



These are the usual direct and exchange integrals which can be written more explicitely as 
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In this way eq. @ becomes 

H = H SSH + - ^2 ^2 Vaa^a'-ja^SaC^C^C-ya/CSa; (6) 
aftjS aa' 

where the prime on the first summation is to omit the terms with a = 5, (3 = 7 and = 7, 
(3 = 8 and 

H SSH = £ an acT + ]r [(Ua/3 - J a p) fi^h^ + U^n^hp^} . (7) 

Off aj3a 

The relevance of this formal partition is twofold: i) it naturally leads to a perturbation ex- 
pansion in the off-diagonal interactions which are in general smaller than the semi-diagonal 
ones; ii) moreover, the unperturbed term 7i SSH is one-body like, with single Slater deter- 
minants as exact eigenstates. The SSH approach defined in Ref. consists in assuming 
TH ~ 7{ SSH , which amounts to neglecting the second and higher order contributions in the 
off-diagonal interactions, the first order one being exactly zero. 

The assumption that the off-diagonal V at p nj s are negligible with respect to the semi- 
diagonal ones is implicit in all the methods which describe electron correlation in terms of 
the Hubbard model, either in its original form,! including only on-site interaction between 
opposite spin electrons — proportional to U a p — , or adding the interaction between parallel 
spin electrons as well, proportional to {U a p — J a p). The important point here is that when the 
Hubbard model is applied to an isolated QD, i.e. to a single site, the Hubbard Hamiltonian 
turns out to be one-particle like: this is so because the inter-site hopping of the traditional 



H SSH n 

1 l , 1 ^ao- 



is zero. As a 



Hubbard Hamiltonian is absent in this case and the commutator 
consequence, the Slater determinants, eigenstates of the single particle Hamiltonian H , are 
exact eigenstates of H SSH as well. 

Within the SSH approach the total energy of N electrons in a QD structure is given by 



E SSH (N) = <V| H SSH \$ N ) =J2 £ a(n a a) + \ £ Papfo-*) + (U a p - J a p) (hp.)] (n aa ) , 

aicr Z after 

(8) 

where <$> N \ is a Slater determinant, eigenvector of Hq, and ( ) denotes the average over 
the many-particle eigenstate, which in our case simply reduces to the orbital occupation 
number. 

The proposed SSH approach shares in common with Hartree-Fock methods the form of 
the total energy, which in both schemes is expressed as the average of the exact Hamiltonian 
over a single Slater determinant; the variational prescription — allowing for the construction 
of optimal single-particle orbitals through the selfconsistent solution of a single-particle 
eigenvalue problem — is not present in the SSH approach. We notice however that the 
importance of self-consistency is strongly related to the relative weight of Coulomb matrix 
elements: the HF potential entering the self-consistent HF one-particle Hamiltonian is in 
fact related to the direct and exchange Coulomb integrals; similarly the SSH approximation 
is exact — without any need of self-consistency — whenever the higher-order contributions 
from the off-diagonal Coulomb matrix elements is negligible. For this reason we expect that 
a lower localization of the confined single-particle states in 3D with respect to 2D, giving 
rise to smaller non-diagonal Coulomb integrals, will reduce the difference between HF and 
SSH results. To check in detail this point we have explicitely performed HF calculations; we 
have used in particular the matrix form of the unrestricted HF equation.! 

Whenever possible, it is obviously useful to compare the outcomes of different approx- 
imate schemes with exact results. This is done in the next section where we consider the 
exactly solvable two-electron QD (artificial Helium) in different confinement regimes; we will 
show that the differences between HF and SSH results will be always comparable with those 
between HF and exact results and that they scale with the dimensionality of the confining 
potential. 
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III. THE TWO-ELECTRON PROBLEM 



In this section we will study the motion of two electrons within a QD structure in two 
and three dimensions. In this case, the exact Hamiltonian reduces to 

H = H Q {1) + H (2) + — ^ 1 , (9) 

k \r t - r 2 \ 

Here is the position of the electron, T*j = (xi,yi) in 2D or r» = (x^y^Zi) in 3D, and p t 
the corresponding momentum. 

To solve this equation, we perform the standard transformation^ to center of mass 
(CM) coordinates, R = (n + r%) /2, P = p l + p 2 , and relative-motion (rm) coordinates, 
r = r\ — r 2 , p = (pi — p 2 ) /2. The two-body Hamiltonian thus splits into a CM and a rm 
part: 

H = H CM + Hrm, (10) 

where 

v 2 1 e 2 

ffrm = f + -^ 2 r 2 + — , (12) 

2/i 2 /tr 

with M = 2m*, and /x = m*/2. The CM hamiltonian -ffcM has the form of a simple harmonic 
oscillator. For the rm Hamiltonian, H rm , it is easy to separate variables and obtain a radial 
differential equation, which gives solutions with the same set of quantum numbers as for 
the harmonic oscillator. Solutions and notations for the 2D and 3D case are summarized in 
Appendix A for both Hqm and H rra . 

By denoting the CM and rm quantum numbers with capital and small letters respectively, 
the eigenvalues for the two-particle system can be written as 

E NM ,nm = hoJ (2N + \M\ + 1) + e nm (13) 

in the 2D case, and 




(14) 



in the 3D spherical case, the cylindrical 3D Helium QD reducing to an effective 2D one (see 
Appendix A). Here e nm and e n e are the rm eigenvalues in 2D and 3D, respectively. Note that 
degeneracy is strongly reduced by Coulomb interaction with respect to the non-interacting 
case. 

The corresponding two-particle total eigenfunctions are 



for a 3D sphere. Here $(-R) and <p(r) are respectively the spatial CM and rm eigenfunctions, 
and x(5, S z ) is the spin function of a state with total spin h 2 S (5 + 1) and z projection 
S z . Note that the parity of the rm spatial eigenfunction is defined (total orbital angular 
momentum and spin are conserved) and connected with the value of total spin by the 
antisymmetry of the two-particle total wavefunction ^(r\, sx',r%, S2). For both the disk and 
the cylinder, this implies that if m is even, the state is a singlet (S = 0), and if m is odd, 
the state is a triplet (5 = 1). Similarly for the sphere case, if I is even, the state is a singlet 
(5 = 0), and if £ is odd, the state is a triplet (5 = 1). 

In the above eigenvalues and eigenfunctions of the two-electron dot, the ingredients re- 
lated to the CM Hamiltonian are known analytically (see Appendix A), while the rm energies 
and wavefunctions must be determined numerically. This is done by exact diagonalization 
of the rm eigenvalue problem (Appendix A), thereby yielding the full 2D and 3D spectrum 
of the QD helium. 

Before comparing these exact results with the SSH approach, we point out that the 
Hamiltonian @ can be translated into a second-quantized form; this is done in terms of the 
same quantum numbers using CM and rm variables. The two-particle Hilbert space is the 



(r) x(S, S z ) ; 



(15) 



for 2D and the 3D cylinder, and 



^NLM s ,Tdm z ;SS t {ri, Si, r 2 , S 2 ) = &NLM Z (R) ^nlm z { T ) X(5, S z ) J 



(16) 



S 



Kronecker product of CM and rm single-particle spaces, generated respectively by the basis 
{|^V")}jv (with eigenvalues E N and creation operators ar N ) and {|^)} n (with eingenvalues e n 
and creation operators a n ). For simplicity, N and n label here the whole set of CM and rm 
quantum numbers, respectively The second-quantized form of the two-particle hamiltonian 
7i, in this variables, is then given by 

H = ^2 E N a} N a N + ^nMi^n + V nn >a) n a nl . (17) 

N n nn ' 

This formulation allows us to obtain the result of the previously discussed Hubbard model, 
by simply neglecting all off-diagonal matrix elements in equation (|17|): 

H SSH = ]T E N a\a N + £ (e„ + V nn ) a n a n , (18) 

TV n 

i.e. manifestly the non interacting Hamiltonian "renormalized" by Coulomb interaction. 

In order to check the reliability of the approximations and the role of dimensionality of 
the confining potential we have calculated ground state properties for QD's with different 
confinement energies, i.e. different values of hu>o, assuming either a 2D or a 3D confining 
potential. The quality of the ground-state eigenfunctions can be probed by the spatial pair 
correlation function f(r): 

f(r) = Kfes(r-r i + r j )). (19) 
Vi / 

Because of the circular symmetry, f(r) depends only on the modulus of the relative distance 
r. Here, the factor K is chosen in such a way that, if we define the dimensionless relative 
distance x = r^j2m*uj /h, the quantity g(x) = xf(x) for the 2D and the 3D cylinder case, 
and g(x) = x 2 f(x) for a 3D sphere is normalized: 

POD 

/ g(x) dx = 1. 
Jo 

We have calculated this quantity both exactly and according to the SSH scheme for a in- 
plane confining energy huo=5 meV (throughout the paper we use m* = 0.065 m e inside 
the dot and m* = 0.079 m e outside, k = 12.98, as in the QD of ref. § m e is the electronic 
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mass); the results are shown in Fig. [I]. The deviations between SSH and exact results clearly 
depend on the dimensionality of the confining potential: in disk-shaped 2D QD's the SSH 
approximation is found to overestimate the probability of finding the two electrons close 
together, in analogy with HF results;! the differences between exact and SSH results are 
significantly reduced assuming a 3D confining potential. This result is coherent with what is 
found for other ground-state properties: Fig. ||(a) shows the ground-state energies calculated 
for dots with different confinement energies, huo, in the range between 4.5 and 10 meV. We 
compare the exact results with the outcomes of HF and SSH calculations assuming 2D and 
3D confinement potentials. Notice that the differences between HF and SSH are always 
smaller — by approximately 50% — than the corresponding differences with respect to the 
exact results; moreover the 3D confinement reduces the overall deviation of both HF and 
SSH by about 60%. 

Since the SSH scheme is exact at the first perturbative order in the off-diagonal ma- 
trix elements of the e-e interaction, it is interesting to check the importance of the next 
perturbative corrections. Details of how the perturbative expansion is actually performed 
for the helium QD are reported in Appendix B. Figure |3] reports exact and SSH ground- 
state energies compared with the results of second order perturbation theory, showing that 
second-order corrections become much smaller if a 3D confinement is assumed. 

The situation becomes more complicated when considering dots with smaller confinement 
energies: in this case the HF and SSH differ from the exact result not only quantitatively 
but also qualitatively, predicting the two-particle ground state to be a triplet instead of 
a singlet, as it should be. This is shown in Fig. 0(b) where again the exact, HF and 
SSH results are shown for dots of different confinement energies. The difference between 
triplet and singlet state energies decreases with increasing confinement energy both for SSH 
and HF approximations until a crossover occurs; assuming a 3D confining potential the 
confinement energy of this crossover is reduced, and this again is true both for SSH and HF 
appr oximat ions . 

We may summarize this analysis on helium QD by concluding that the assumed 3D con- 

10 



finement potential reduces the differences between approximate (SSH and HF) solutions with 
respect to the exact ones, both in terms of ground-state eigenfunctions and eigenenergies. 



IV. THE MANY-ELECTRON PROBLEM 

The key quantity that characterizes single-electron transport into a QD is the addition 
energy, i.e., the energy A(N) required to place an extra electron into a dot that is initially 
occupied by iV — 1 electrons. Such quantity, analogous to electron affinity in atomic physics, 
can be measured experimentally as a function of N. It has been shown! that the measured 
voltage increment AA between successive single-electron tunneling processes — i.e., between 
two successive maxima in the conductance — peaks at "magic" values of N corresponding 
to the filling of complete shells (N = 2,6, 12), as well as to half-shell filling (e.g. N = 4). 
The existence of these half-shell filling features is reminiscent of Hund's rule in atomic 
physicsJII! and is intimately related to electron-electron interaction. 

The results of SSH theory for the addition-energy variations, AA(N) = A(N + 1)—A(N), 
are displayed in Fig. |] as a function of the electron number N for two different 3D cylindrical 
quantum dots. Here, A(N) is obtained as E SSH (N) - E SSH (N - 1), where E SSH (N) is 
the ground-state energy in Eq. As we can see, AA(N) exhibits peaks corresponding 
both to complete and half shell filling, thus well reproducing the experimental evidence in 
Ref. 0. This behaviour is the result of the interplay between single-particle contributions 
and electron-electron repulsion: the single-particle term favors complete shell filling, while 
the repulsion among parallel-spin electrons, smaller than the repulsion among opposite-spin 
ones, makes the configurations with maximum total spin energetically favored (Hund's rule). 
This is the physical origin of the half shell-filling structure: indeed, adding an electron to 
a half-filled shell forces the double occupancy of a level; consequently, AA is raised by the 
dominant Coulomb repulsion U aa between opposite-spin electrons on the same level. 

For some non-closed shell configurations the total spin turns out to be not determined 
by Hund's rule: in particular, for N=1Q we find a ground state with total spin S=0. Sim- 
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ilar deviations from Hund's rule have been found for large electron numbers (N >20) and 
associated to spin-density-wave instabilities; for smaller numbers (iV=16,18) the same 5=0 
spin-density- wave state has been found to be a low-energy "spin isomer" , slightly higher in 
energy than the ground state configuration.^ 

From our calculations we may say that the S=0 configuration in dots with large N may 
be favoured by the reduced repulsion between electrons in high shells: in the fourth shell, 
for instance, the Coulomb integrals U a p relative to orbitals with higher values of the orbital 
momentum may be smaller than the corresponding terms relating two levels with smaller 
angular momentum; the double occupation of an orbital with high orbital momentum m 
(i.e., the level with n=0, m=3, see App. A for the notation) with antiparallel-spin electrons 
may therefore cost less than having parallel-spin electrons on different degenerate orbitals, 
but with smaller m (i.e., the level n=l, m=l). The same interplay also explains the peaks 
in AA(N) for iV =14 and iV=18.0 

We want to stress that also in the case of many electrons the reliability of the results of 
SSH approach is comparable with HF ones. The explicit comparison between the addition 
energy variation calculated according to SSH and HF schemes and for 2D and 3D confine- 
ments is reported in Fig. [|, showing that AA always peaks at the same electron numbers 
and that the agreement between SSH and HF results improves on going from the 2D to the 
3D confinement model. 

Ground-state configurations and filling rules change when a magnetic field is applied. It 
affects both single particle energies and Coulomb and exchange integrals through the induced 
changes in the wavefunction localization. Figure |6] shows the U- and J- integrals vs. B for 
the first states, obtained for huo=7.5 meV. For comparison, we also show the corresponding 
quantities calculated within a strictly 2D-confinement model. We can see that {/-integrals 
describing the interaction between opposite-spin electrons are a few meV smaller in the case 
of 3D confinement, while the differences in the interaction between parallel-spin ones are 
much smaller. This is going to affect dramatically the energy balance which determines 
ground-state configurations, thus clearly showing the failure of a pure 2D description of 
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state-of-the-art QD structures. 

As already mentioned, according to the SSH approach the off-diagonal matrix elements of 
the electron-electron interaction are assumed to be negligible. In Fig. || the values of two of 
them are reported as functions of the applied magnetic field. As expected, we clearly see that 
for any B value they are negligible compared to all the other semi-diagonal contributions, 
and even more so in 3D with respect to the 2D case. 

Figure shows the total energy E as a function of the applied magnetic field B for 
different values of the electron number N in a dot with confinement energy huo=7.5 meV. It 
appears that for sufficiently large values of B the Hamiltonian term linear in the magnetic 
field becomes dominant, making configurations with higher total angular quantum number 
energetically favorable. This is also the physical origin of the wiggles in the A(N) vs. B plot 
shown in Fig. |8] and observed in the experiments reported in Ref. |2|. 

Other authors have explicitely considered the question of dimensionality in theoretical 
modelization of semiconductor QDs. Kumar et aM computed self-consistently the one- 
particle confining potential in a square QD. According to their results, our assumption of an 
in-plane parabolic confining potential plus a well in the perpendicular direction is seen to be 
quite reliable and general, as well as the ansatz of considering only the ground state motion 
along z, at least for few electron dots. Steinebach et aZ.0 pointed out the importance of 
a full 3D model to treat spin density excitations (SDEs) in semiconductor wires and dots. 
Specifically, they used the analogous of eq. (|A4[) as effective 3D Coulomb interaction, and 
they found that a 2D description artificially enhances the interaction strenght and is unable 
to predict experimental Raman spectra. The necessity of a 3D modelization is then seen to 
emerge not only in the description of ground state and single particle processes, like addition 
spectra, but also in two particle processes, like SDEs. 
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V. SUMMARY AND CONCLUSIONS 



We have presented a theoretical investigation of Coulomb-correlation effects in semicon- 
ductor quantum dots. In particular, we have performed a detailed analysis of the addition- 
spectrum problem for few-electron quantum-dot structures (macroatoms), pointing out pos- 
sible analogies with more conventional Coulomb-correlation effects in atomic physics. 

Our primary goal was to understand to which extent the various approximation schemes, 
such as Hartree-Fock or Hubbard models, are able to properly describe Coulomb correlation 
in realistic, state-of-the-art QDs. To this end, we have first compared approximate results to 
the exact solution for the prototypical case of a two-electron system, the so-called quantum- 
dot Helium; we have repeated such analysis for different dimensionalities, considering the 
3D (spherical and cylindrical geometry) and the pure 2D structures. The main result is that 
the degree of accuracy of any approximation scheme depends strongly on the dimensionality 
of the problem. More specifically, the pure 2D model — often used for the description 
of quantum dots — is found to give approximated results which differ significantly from the 
exact solution. We have demonstrated that this is not a general failure of the approximation 
scheme, but it reflects a rather pathological behaviour originating from the unphysical nature 
of the pure 2D model. Indeed, for the case of a 3D cylindrical model — which provides a 
much better description of realistic QD structures — the difference between exact solution 
and approximated results is found to be much smaller, thus confirming the validity of the 
various approximation schemes considered. 

The same analysis has been then extended to many-electron systems for which addition- 
spectra measurements are available. Using different approximation schemes, we find that 
the deviations between the full 3D and the simplified 2D quantum-dot model are very 
significant. The full 3D model is found to reproduce the experimental data for a large 
class of QD structures where simplified 2D models fail. We conclude that this is due to 
the unphysical character of the pure 2D confinement for which the various approximation 
schemes often yield unreliable results. A proper description of the QD structure in terms 
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of fully 3D single-particle wavefunctions is therefore required; we have shown that in this 
case approximate approaches can give an accurate description of correlation effects in the 
macroatoms made available by present semiconductor technology. 

We are grateful to D. Pfannkuche and D. Vanossi for useful discussions. This work was 
supported in part by the EC through the TMR Network "Ultrafast quantum optoelectron- 
ics". 



APPENDIX A: CM AND RM SOLUTIONS FOR THE TWO-ELECTRON DOT 



The present appendix is organized as follows: In Sec. |A 1| we shall show how to reduce 



the 3D-cylindrical Helium problem to an effective two-dimensional one; In Sect. [A2| we shall 
summarize the 2D and 3D solutions of the one-particle Schrodinger equations for the center- 
of-mass and the relative-motion Hamiltonians, i^cM and H T m, as defined in eqs. (|TT| ) and 

(0)- 



1. 3D eigenvalue equation for the cylindrical QD 

If only the lowest single-particle state 4>o{z) of the quantum well is relevant to the two- 
electron motion, we can write the spatial-part ^(Vi; r 2 ) of our Helium wavefunction as 

r 2 ) = ^{xi,yi; x 2 , y 2 ) 0oOi) </>o{z 2 ) • (Al) 

This approximation is well justified for most cases of interest. Indeed, for the typical QD 
structure used in the experimental investigation of addition spectral (quantum-well width 
L = 12 nm, barrier height Vq = 200 meV, the energy separation between ground and first 
excited state along z is 56 meV, about one order of magnitude larger than typical in-plane 
single particle confinement energies. 

Let us now consider the global Schrodinger Equation corresponding to the exact Helium 
Hamiltonian of Eq. (Q) 

H^(r 1 ;r 2 ) = E^(r 1 ;r 2 ); (A2) 
15 



by substituting eq. multiplying both sides by 4>q(zi) ^0(^2) and integrating over z\ and 

z 2 , we obtain: 



e 2 



2e z + J2H (i) + -c(\r 1 -r 2 \) 



^(r i; r 2 ) =^(r i; r 2 ). (A3) 



i=l 

The eigenvalue equation is then reduced to a 2D one, since = (xj,t/j) and c(r) = 
c (\ r i — r 2\) is an effective Coulomb potential, accounting for the geometry of the system: 

c(r) = r^jj^dm^L; ( A4) 

J -00 J -00 s J r 2 + ( Zl _ Z2 ^ 

From now on we will drop the constant ground-state energy along z (eg). As a first step, we 
evaluate c(r) by solving the quantum well eigenvalue problem (allowing for different values 
of the effective mass in the well and in the barrier); Then, we numerically integrate eq. ( |A 1| ) . 
It is easy to show analytically some important properties of c(r), namely that 

< rc(r) < 1 Vr, (A5) 

lim rc(r) = 0, (A6) 
) — >o 

lim rc(r) = 1. (A7) 

r— *oo 

These properties tell us that for large distances r c(r) tends to the bare Coulomb potential, 
and that it is however strongly reduced in the neighborhood of the origin, i.e., the more 
relevant space region in the computation of Coulomb and exchange integrals. Figure ^| shows 
such effective Coulomb potential C multiplied by the dimensionless variable x (introduced 
below) as a function of x for different values of the quantum-well width: A monotonous 
behavior is apparent, going from the bare Coulomb-potential value in the zero- width limit 
(the function is constant and equal to 1), into progressively smaller values, towards the 
infinite-width case. 
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2. Exact solutions 



Let us first consider the CM equation, which has the form of a standard harmonic oscil- 
lator and can thus be solved analytically. For the 2D case (3D cylindrical), its eigenvalues 



arc: 



e™ M = hu (2N+\M\ + 1) 
N = 0,1,2,... M = 0,±l,±2,... (A8) 



and the corresponding ortonormalized eigenfunctions (the so called "Fock-Darwin" statesu't3) 



arc 



*&r„(r,s) = (s,r\NMa) = A^ - ( ^ |M|) , e^r^^L^ (Ar 2 ) X(T (s) . (A9) 

In the 3D spherical case,@ the eigenvalues are 

ef L = hu (2N + L + 
N = 0,1,2,... L = 0,1,2,... (A10) 

and the ortonormalized eigenfunctions ^NLM z A r i s ) are 



$f ML Jr,s) = (s,r\NLM z a) 



2XL+lNl r^L^{^)Y LM ^)x,{s). 



\ t(n + l + 

(All) 

Here, A = m *^° , L P N are generalized Laguerre polynomials,^ Y is the usual Gamma function, 
Xo- denotes the spin function, and Y LMz are the spherical harmonics. We have used polar 
coordinates throughout: r = (r, </?) in 2D (3D cylindrical), case and r = (r,d,(p) in the 3D 
spherical case. For the 2D (3D cylindrical) the quantum numbers are (N,M,a): N is the 
radial quantum number, M the angular momentum quantum number (in this case the total 
angular momentum coincides with the component along z, L z = —hM), and a the spin 
component along z. In the 3D spherical case, on the other hand, the quantum numbers are 
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given by (N, L, M z , a): here L is the total angular momentum quantum number, and M z is 



the magnetic quantum number, M z = —L, —L + 1, . . . , L. 

Let us now come to the single-particle Schrodinger equation for the rm Hamiltonian of 
Eq. flL^D- In this equation the variables are easily separable, and the problem is reduced 
to the solution of a radial differential equation. For the 2D (3D cylindrical) case, the rm 
eigenfunction in coordinate space is 



-mup 



if, 



Rn.m.{r) 



^7T 



(A12) 



where R nm (r) is the solution of the radial Schrodinger equation 



d 2 R nm (r) . ldRnmir) 



dr 2 ^ r dr 

we have employed the notations A = fiuo/h, a = 2/ie 2 //t/i 2 , and k 2 m = 2fie nm /h 2 , where e nm 
is the rm eigenvalue. The effective Coulomb potential c(r) is simply 1/r in the 2D case and 
it is defined in Sect. [A 1| for the 3D cylindrical case. 

For the 3D spherical case, the rm eigenfunction in coordinate space is 



+ 



Z..2 

« ... 



m 



X r — ac[r) — 



0; 



(A13) 



(A14) 



with R n i{r) satisfying the radial eigenvalue equation 

}2 D j„\ oad .u\ r a £(£ + 1) 



O-R.nir) | 2 dR ne (r) | 



RnAr) 



0; 



(A15) 



where, again, we put k 2 e = 2fie n £/h 2 , and e„£ is the rm eigenvalue. 

In order to obtain an exact solution for the rm eigenvalue problems, we rewrite Eqs. (|A13|) 
and (|A15|) in terms of the dimensionless variable x = A 1 ' 2 ?". For the 2D (3D cylindrical) 
case, Eq. ( |A13j ) becomes 



_d_ f dR nm (x) \ 
dx \ dx I 



in 



(xx C (x) ~\~ k nm x x 

X 



Rum (%) 



0. 



Rnmy^l RumiX) 
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C(x) = \- 1/2 c(\- 1/2 x) (A16) 
(for the 2D case it is simply C{x) = 1/x), while for the 3D case Eq. ( |A15| ) transforms into 



dx 2 



£(£ + 1) a 



+ k„n X 



X 2 X 



Xni(x) = 0, 



Xne(x) = Xm(r) Xne(r) = Rni T . (A17) 



The dimensionless parameters are at = X^^a = 2ylZ* /huj , k 2 = k 2 /X = 2e a /huj; 1Z* = 
e 4 m* /2n 2 h 2 is the effective Rydberg energy. Actually, exact analytic solutions exist, but they 
are limited to 2D and 3D spherical cases only|iB thus, we have chosen to solve Eqs. (|A16|) 
and ( A17|) by standard numerical methods. We stress that the numerical accuracy depends 



on the accurate specification of the boundary conditions, that we impose through analytical 
asymptotic formulae for eigenfunctions near to the singular points and +oo, following the 
general methods of Ref. [22|. In this way the numerical solution is very stable and efficient, 
thus overcoming possible difficulties related to the singlet ground state;! in out calculations 
energy values are obtained with a nominal relative error of the order of 10 -8 . 

APPENDIX B: HELIUM PERTURBATION THEORY 

We employ the standard Rayleigh-Schrodinger perturbation theory to correct the SSH 
eigenvalues in Eq. (|T8|) at the second order in the off-diagonal Coulomb matrix elements 
entering the total Hamiltonian flTTp. In the remaining part of this section we shall consider 
the 2D and 3D spherical cases, by neglecting the center-of-mass motion. 

For the 2D case, the rm SSH eigenvalues e££ H are given by0 



= Hujo (2n + \m\ + 1) + y/K*hu) S 2D (n } m) , (Bl) 

S*> (n , m) "H + "f "'(-l>-"«2. + !)!!] VI' \. (B2) 

' ^ {n-s-l)\2 2s + 2 [{s + l)\] 2 {\m\ + s + l)\) K ' 

while for the 3D case we have 
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rNm 






m 


I 



e s J H = nu (2n + e+^j + ^Jn*huj S m (n, £) , (B3) 

2^r 2 (£+|) \ s tS2^+3(n- S -l)!r( S +l+|)[(s + l)!] 2 J 
The first-order correction due to non-diagonal Coulomb matrix elements is equal to zero. 

(2) 

The second-order correction A^q to the ground-state energy is given by the well-known 
expression 

A£ 2) = S F SSH _ ZSSH ( B5 ) 

(see notation in eq. (|17D ). The idea now is to look for analytic expressions for the off-diagonal 
integrals Vo n and then to perform a numerical summation. However, expressions like those 
obtained in Eqs. (|B2| )- (|B4| )ii are not useful, since each integral is given by an alternated-sign 
summation and numerical errors become rapidly critical as the quantum number n increases. 
In contrast, the solution can be obtained using an integration trick suggested in Ref. 24), so 
that all the terms in the summation are obtained with the same sign. For the 2D case one 
gets 

r(n + |) 



Vno = \ lZ*hca, 



1 



r(n+l) 



V n0 w \IK*hujQ — n -> oo. (B6) 

W2 



For the 3D case one obtains 



V n0 = sJ1Z*TujJq 







TiY(n + 1) T 


[n + 


1! 


r 





„n ~ \/TZ*tuv i — n -»• oo. (B7) 

As already pointed out, now the generic terms (|B6p-(|B~7|) in the sum ( |B5| ) have the same 
sign and the summation can be easily performed. The result for the 2D case is 
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A4 ' = -TV (0.691), (B8) 

and for the 3D case 

A4 2) = -TV (0.156). (B9) 

Note that the 3D term is significantly smaller than the corresponding 2D one, and that in 
the 3D case the series converges faster. 
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FIGURES 

FIG. 1. Ground-state spatial pair-correlation function g(x) for the 3D (spherical) and 2D 
two-electron QD: exact and SSH results are reported. Here, x = ^ 2m £° r ^ 2 is the dimension- 
less relative radial coordinate and g{x) is normalized in such a way that f£°g(x) dx = 1. The 
in-plane confinement energy is %ojq = 5 meV. 

FIG. 2. (a) Ground-state energy of the artificial He QD as a function of the confinement energy 
huo, calculated within different approaches (exact, SSH, and HF). The range of Tiloq is 4.5-^10 meV. 
The two panels correspond to 2D and 3D (cylindrical) geometries. The ground-state configuration is 
always a spin-singlet, (b) Spin-singlet (spin-triplet) energies versus confinement energies Tuvq. The 
range of Tlujq is l-f-4.5 meV. The two panels are relative to 2D and 3D (cylindrical) geometries. The 
exact ground state is always a singlet, while a singlet-triplet crossover occurs for both approximated 
schemes in the low-energy region. 

FIG. 3. Ground-state energy of the artificial He QD as a function of the confinement energy 
Tiuo, as obtained via exact diagonalization, SSH approximation scheme, and Rayleigh-Schrodinger 
perturbation theory at the second order in the off-diagonal Coulomb matrix elements. Both 3D 
(spherical) and 2D cases are shown. 

FIG. 4. Calculated SSH addition-energy increment /S.A as a function of the total number N 
of electrons for two different QD structures, both characterized by a parabolic potential in the xy 
plane (confining energy fuvo) and by a finite-barrier quantum-well potential along the z direction 
(3D cylindrical model). 

FIG. 5. Comparison between SSH and HF addition-energy increments AA as a function of the 
total number N of electrons in the dot. Here, the upper panel corresponds to the 2D geometry 
while the lower one corresponds to the 3D cylindrical model. The in-plane confinement energy is 
hojo = 7.5 meV. 
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FIG. 6. Coulomb (U a -p) and exchange {J a -p) integrals as well as off-diagonal Coulomb matrix 
elements {V a ^- 1: s), as functions of the magnetic field B for both 2D and 3D cases. Here, a and (5 
denote the sets of radial and angular quantum numbers (n, m) for the various single-particle states 
involved in the two-body interaction process. The in-plane confinement energy is hcoo = 7.5 meV. 

FIG. 7. Total energy E SSH as a function of the applied magnetic field B corresponding to 
N electrons in a dot with confinement energy Tiuiq = 7.5 meV. For any given value of N, all the 
possible configurations, denoted by the usual atomic physics terms 2S+l L, have been considered. 

FIG. 8. Addition energy A(N) as a function of the magnetic field B calculated for a realistic 
(3D) QD structure with confinement energy Tiujq = 7.5 meV and for different values of N. The 
labels indicate the electronic terms for the ground-state configurations, that depend on B. 

FIG. 9. Plot of the effective Coulomb potential C(x) multiplied by the dimensionless coordi- 
nate x (= xC(x)) as a function of x for different values of the quantum well width L and for a 
confinement energy Tiujq = 5meV. Notice that in the limit L — ► (2D case) C(x) — ► \jx and, 
therefore, xC(x) — > 1. 
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